﻿<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META content="text/html; charset=utf-8" http-equiv=Content-Type></HEAD>
<BODY><PRE>double addQuad(double a, double b, double c=0.0, double d=0.0)
{
	return sqrt(a*a+b*b+c*c+d*d);
}

void XsecSumPlotFull(int log = 1)
{
	
	// Data
	//
	
	enum modes{
		kJet,
		kHF,
		kW,
		kZ,
		kWg,
		kZg,
		kWW,
		kttbar,
		kWZ,
		kST,
		kZZ,
		kHWW,
		kNP,
		kNModes,
	};

	int firstBin=kW;
	
	double cdfvals    [kNModes][3];
	double cdfvals_sm [kNModes][3];
	int    cdfisprelim[kNModes];
	double d0vals     [kNModes][3];
	double d0vals_sm  [kNModes][3];
	int    d0isprelim [kNModes];
	
	// Jet
	cdfisprelim[kJet] = 0;
	cdfvals    [kJet][0] = 1.0*1000*1000*1000; 
	cdfvals   [kJet][1] =  0.001; 
	cdfvals   [kJet][2] = 0.;
	cdfvals_sm [kJet][0] = 1.0*1000*1000*1000; 
	cdfvals_sm[kJet][1] =  0.;    
	cdfvals_sm[kJet][2] = 0.;
	d0isprelim [kJet] = 0;
	d0vals     [kJet][0] = 1.0*1000*1000*1000; 
	d0vals    [kJet][1] =  0.001; d0vals    [kJet][2] = 0.;
	d0vals_sm  [kJet][0] = 1.0*1000*1000*1000; d0vals_sm [kJet][1] =  0.;    d0vals_sm [kJet][2] = 0.;
	
	cdfisprelim[kHF] = 0;
	cdfvals    [kHF][0] = 0.1*1000*1000*1000; cdfvals   [kHF][1] =  0.001; cdfvals   [kHF][2] = 0.;
	cdfvals_sm [kHF][0] = 0.1*1000*1000*1000; cdfvals_sm[kHF][1] =  0.;    cdfvals_sm[kHF][2] = 0.;
	d0isprelim [kHF] = 0;
	d0vals     [kHF][0] = 0.1*1000*1000*1000; d0vals    [kHF][1] =  0.001; d0vals    [kHF][2] = 0.;
	d0vals_sm  [kHF][0] = 0.1*1000*1000*1000; d0vals_sm [kHF][1] =  0.;    d0vals_sm [kHF][2] = 0.;
	
	cdfisprelim[kW]    = 0;
	cdfvals    [kW][0] = 2.749*1000; 
	cdfvals   [kW][1] = addQuad (0.01,0.053,0.165)*1000; cdfvals   [kW][2] = addQuad(0.01,0.053,0.165) *1000;
	cdfvals_sm [kW][0] = 2687;       cdfvals_sm[kW][1] =  54;     cdfvals_sm[kW][2] = 54;
	d0isprelim [kW]    = 0;
	d0vals     [kW][0] = 2865.2;     d0vals    [kW][1] = addQuad (8.3,62.8,40.4,186.2);   d0vals    [kW][2] = addQuad (8.3,62.8,40.4,186.2);
	d0vals_sm  [kW][0] = 2687;       d0vals_sm [kW][1] =  54;                             d0vals_sm [kW][2] = 54;
	
	cdfisprelim[kZ] = 0;
	cdfvals    [kZ][0] = 254.9; cdfvals   [kZ][1] = addQuad (3.3,4.6,15.2); cdfvals   [kZ][2] = addQuad(3.3,4.6,15.2);
	cdfvals_sm [kZ][0] = 251.3; cdfvals_sm[kZ][1] =  5.0;                   cdfvals_sm[kZ][2] = 5.0;
	d0isprelim [kZ] = 1;
	d0vals     [kZ][0] = 264.9; d0vals    [kZ][1] = addQuad (3.9,8.5,5.1,17.2); d0vals[kZ][2]    = addQuad(3.9,8.5,5.1,17.2);
	d0vals_sm  [kZ][0] = 251.3; d0vals_sm [kZ][1] =  5.0;                       d0vals_sm[kZ][2] = 5.0;
	
	cdfisprelim[kWg] = 1;
	//	cdfvals    [kWg][0] = 19.11; cdfvals   [kWg][1] = addQuad (1.04,2.40,1.11); cdfvals   [kWg][2] = addQuad(1.04,2.40,1.11);
	cdfvals    [kWg][0] = 18.03; cdfvals   [kWg][1] = addQuad (0.65,2.55,1.05); cdfvals   [kWg][2] = addQuad(0.65,2.55,1.05);
	cdfvals_sm [kWg][0] = 19.3 ; cdfvals_sm[kWg][1] =  1.4;                     cdfvals_sm[kWg][2] = 1.4;
	d0isprelim [kWg] = 0;
	d0vals     [kWg][0] = 14.8; d0vals   [kWg][1] = 2.1; d0vals   [kWg] [2] = 2.1;
	d0vals_sm  [kWg][0] = 16.0; d0vals_sm[kWg][1] = 0.4; d0vals_sm[kWg] [2] = 0.4;
	
	cdfisprelim[kZg] = 1;
	//	cdfvals    [kZg][0] = 4.9; cdfvals   [kZg][1] = addQuad (0.3,0.3,0.3); cdfvals   [kZg][2] = addQuad(0.3,0.3,0.3);
	//	cdfvals_sm [kZg][0] = 4.7; cdfvals_sm[kZg][1] =  0.4;                  cdfvals_sm[kZg][2] = 0.4;
	cdfvals    [kZg][0] = 4.6; cdfvals   [kZg][1] = addQuad (0.2,0.3,0.3); cdfvals   [kZg][2] = addQuad(0.2,0.3,0.3);
	cdfvals_sm [kZg][0] = 4.5; cdfvals_sm[kZg][1] =  0.3;                  cdfvals_sm[kZg][2] = 0.3;
	d0isprelim [kZg] = 0;
	d0vals     [kZg][0] = 4.96; d0vals   [kZg][1] =  addQuad(0.3,0.3);  d0vals   [kZg][2] = addQuad(0.3,0.3);
	d0vals_sm  [kZg][0] = 4.74; d0vals_sm[kZg][1] = 0.22;               d0vals_sm[kZg][2] = 0.22;
	
	cdfisprelim[kWW] = 1;
	//	cdfvals    [kWW][0] = 13.69; cdfvals   [kWW][1] = addQuad (2.32,1.56,0.119); cdfvals   [kWW][2] = addQuad(2.32,1.56,0.119);
	cdfvals    [kWW][0] = 12.1; cdfvals   [kWW][1] = addQuad (0.9,1.6); cdfvals   [kWW][2] = addQuad(0.9,1.4);
	cdfvals_sm [kWW][0] = 11.66;  cdfvals_sm[kWW][1] =  0.7;                      cdfvals_sm[kWW][2] = 0.7;
	d0isprelim [kWW] = 0;
	d0vals     [kWW][0] = 13.8; d0vals   [kWW][1] = addQuad (4.8,1.2,0.9); d0vals   [kWW][2] = addQuad(3.8,0.9,0.9);
	d0vals_sm  [kWW][0] = 12.4; d0vals_sm[kWW][1] =  0.8;                  d0vals_sm[kWW][2] = 0.8;
	
	cdfisprelim[kttbar] = 1;
	//	cdfvals    [kttbar][0] = 8.2;  cdfvals   [kttbar][1] = addQuad (0.5,0.8,0.5); cdfvals   [kttbar][2] = addQuad(0.5,0.8,0.5);
	cdfvals    [kttbar][0] = 7.0;  cdfvals   [kttbar][1] = addQuad (0.3,0.4,0.4); cdfvals   [kttbar][2] = addQuad(0.3,0.4,0.4);
	cdfvals_sm [kttbar][0] = 6.6 ; cdfvals_sm[kttbar][1] =  0.8;                  cdfvals_sm[kttbar][2] = 0.8;
	d0isprelim [kttbar] = 1;
	d0vals     [kttbar][0] = 8.3; d0vals   [kttbar][1] = addQuad (0.6,0.9,0.5); d0vals   [kttbar][2] = addQuad(0.5,0.9,0.5);
	d0vals_sm  [kttbar][0] = 6.6; d0vals_sm[kttbar][1] =  0.8;                  d0vals_sm[kttbar][2] = 0.8;
	
	cdfisprelim[kWZ] = 1;
	//	cdfvals    [kWZ][0] = 4.3;  cdfvals   [kWZ][1] = 1.1;  cdfvals    [kWZ][2] = 1.4;
	cdfvals    [kWZ][0] = 4.3;  cdfvals   [kWZ][1] = addQuad (1.3,0.4);  cdfvals    [kWZ][2] = addQuad (1.1,0.4);;
	cdfvals_sm [kWZ][0] = 3.65; cdfvals_sm[kWZ][1] = 0.26; cdfvals_sm [kWZ][2] = 0.26;
	d0isprelim [kWZ] = 1;
	d0vals     [kWZ][0] = 2.7;  d0vals    [kWZ][1] = 1.7;  d0vals     [kWZ][2] = 1.3;
	d0vals_sm  [kWZ][0] = 3.65; d0vals_sm [kWZ][1] = 0.26; d0vals_sm  [kWZ][2] = 0.26;
	
	cdfisprelim[kST] = 1;
	//	cdfvals    [kST][0] = 3.0;  cdfvals   [kST][1] = 1.2;  cdfvals    [kST][2] = 1.1;
	cdfvals    [kST][0] = 2.3;  cdfvals   [kST][1] = 0.6;  cdfvals    [kST][2] = 0.5;
	cdfvals_sm [kST][0] = 2.9;  cdfvals_sm[kST][1] = 0.27; cdfvals_sm [kST][2] = 0.27;
	d0isprelim[kST] = 0;
	d0vals     [kST][0] = 4.7;  d0vals    [kST][1] = 1.3;  d0vals     [kST][2] = 1.3;
	d0vals_sm  [kST][0] = 2.9;  d0vals_sm [kST][1] = 0.27; d0vals_sm  [kST][2] = 0.27;
	
	cdfisprelim[kZZ]   = 1;
	//	cdfvals    [kZZ][0] = 0.75; cdfvals   [kZZ][1] = 0.71; cdfvals    [kZZ][2] = 0.54;
	cdfvals    [kZZ][0] = 1.4; cdfvals   [kZZ][1] = 0.7; cdfvals    [kZZ][2] = 0.6;
	cdfvals_sm [kZZ][0] = 1.4;  cdfvals_sm[kZZ][1] = 0.1;  cdfvals_sm [kZZ][2] = 0.1;
	d0isprelim [kZZ] = 1;
	d0vals     [kZZ][0] = 4.3; d0vals   [kZZ][1] = 0.0; d0vals   [kZZ] [2] = 4.3;
	d0vals_sm  [kZZ][0] = 1.6; d0vals_sm[kZZ][1] = 0.1; d0vals_sm[kZZ] [2] = 0.1;
	
	cdfisprelim[kHWW] = 1;
	//	cdfvals    [kHWW][0] = 2.0*0.041/0.1; cdfvals   [kHWW][1] =  0;      cdfvals   [kHWW][2] = 2.0*0.041/0.1;
	//	cdfvals_sm [kHWW][0] = 0.388;         cdfvals_sm[kHWW][1] =  0.0388; cdfvals_sm[kHWW][2] = 0.0388;

	// data from http://www-cdf.fnal.gov/physics/new/hdg/results/hwwmenn_090227/cdf9764_pubnoteHWW36fb.ps
	cdfvals    [kHWW][0] = 1.37*0.439; cdfvals   [kHWW][1] =  0;      cdfvals   [kHWW][2] = 1.37*0.439;
	cdfvals_sm [kHWW][0] = 0.439;         cdfvals_sm[kHWW][1] =  0.0388; cdfvals_sm[kHWW][2] = 0.0388; // errors unknown....
	d0isprelim [kHWW] = 1;
	d0vals     [kHWW][0] = 1.4;    d0vals   [kHWW][1] = 0.0;    d0vals [kHWW][2]    = 1.4;
	d0vals_sm  [kHWW][0] = 0.388;  d0vals_sm[kHWW][1] = 0.0388;  d0vals_sm[kHWW][2] = 0.0388;
	
	cdfisprelim[kNP] = 1;
	cdfvals    [kNP][0] = 0.0001; cdfvals   [kNP][1] = 0;   cdfvals    [kNP][2] = 0.;
	cdfvals_sm [kNP][0] = 0.5;    cdfvals_sm[kNP][1] = 9.5; cdfvals_sm [kNP][2] = 0.4;
	d0isprelim[kNP] = 1;
	d0vals    [kNP][0] = 0.0001; d0vals   [kNP][1] = 0;   d0vals    [kNP][2] = 0.;
	d0vals_sm [kNP][0] = 0.5;    d0vals_sm[kNP][1] = 9.5; d0vals_sm [kNP][2] = 0.4;
	
	
	// Setup
	//
	gStyle-&gt;SetOptStat(0000);
	gStyle-&gt;SetGridColor(15);
	gStyle-&gt;SetPadGridX(1);
	TCanvas* c = new TCanvas("test", "test", 0,0, 600, 600);
	c-&gt;GetFrame()-&gt;SetFillColor(21);
	c-&gt;GetFrame()-&gt;SetBorderSize(12);
	c-&gt;SetBottomMargin(0.10);
	c-&gt;SetRightMargin(0.05);
	c-&gt;SetLeftMargin(0.15);
	c-&gt;SetLogy(log);
	
	// create a 2-d histogram to define the range
	//
	TString titleStr; double minY; double maxY;
	titleStr = "Diboson Cross-Section Measurements at the Tevatron";  minY = 0.1; maxY = 1e5;
	TH2F *h = new TH2F("h", titleStr.Data(), kNModes-firstBin, firstBin, kNModes, 6,  minY, maxY);
	h-&gt;GetYaxis()-&gt;SetTitleSize  (0.04);
	h-&gt;GetXaxis()-&gt;SetTitleSize  (0.04);
	h-&gt;GetYaxis()-&gt;SetTitleOffset(1.45);
	h-&gt;GetYaxis()-&gt;SetLabelSize  (0.04);
	h-&gt;GetXaxis()-&gt;SetLabelSize  (0.04);
	h-&gt;GetXaxis()-&gt;SetLabelOffset(0.0075);
	h-&gt;GetXaxis()-&gt;SetLabelFont  (62);
	h-&gt;SetTitle ("");
	h-&gt;SetXTitle("");
	h-&gt;SetYTitle("Production Cross Section [pb]");
	//  h-&gt;SetTitleOffset(1.2, "Y");
	
	h-&gt;GetXaxis()-&gt;SetBinLabel(1-firstBin+kJet,  "Jets");
	h-&gt;GetXaxis()-&gt;SetBinLabel(1-firstBin+kHF,   "Heavy Flavor");
	h-&gt;GetXaxis()-&gt;SetBinLabel(1-firstBin+kW,    "W");
	h-&gt;GetXaxis()-&gt;SetBinLabel(1-firstBin+kZ,    "Z");
	h-&gt;GetXaxis()-&gt;SetBinLabel(1-firstBin+kWg,   "W#gamma");
	h-&gt;GetXaxis()-&gt;SetBinLabel(1-firstBin+kZg,   "Z#gamma");
	h-&gt;GetXaxis()-&gt;SetBinLabel(1-firstBin+kWW,   "WW");
	h-&gt;GetXaxis()-&gt;SetBinLabel(1-firstBin+kttbar,"t#bar{t}");
	h-&gt;GetXaxis()-&gt;SetBinLabel(1-firstBin+kWZ,   "WZ");
	h-&gt;GetXaxis()-&gt;SetBinLabel(1-firstBin+kST,   "t");
	h-&gt;GetXaxis()-&gt;SetBinLabel(1-firstBin+kZZ,   "ZZ");
	h-&gt;GetXaxis()-&gt;SetBinLabel(1-firstBin+kHWW,  "H#rightarrowWW");
	h-&gt;GetXaxis()-&gt;SetBinLabel(1-firstBin+kNP, "NP?");
	
	h-&gt;Draw();
	
	TAxis* axis2 = h-&gt;GetXaxis()-&gt;Clone();
	axis2-&gt;SetLabelSize  (0.02);
	axis2-&gt;SetLabelOffset(0.01);
	axis2-&gt;SetBinLabel   (1-firstBin+kJet,"");
	axis2-&gt;SetBinLabel   (1-firstBin+kHF,"");
	axis2-&gt;SetBinLabel   (1-firstBin+kW,"");
	axis2-&gt;SetBinLabel   (1-firstBin+kZ,"");
	axis2-&gt;SetBinLabel   (1-firstBin+kWg,"");
	axis2-&gt;SetBinLabel   (1-firstBin+kZg,"");
	axis2-&gt;SetBinLabel   (1-firstBin+kWW,"");
	axis2-&gt;SetBinLabel   (1-firstBin+kWZ,"");
	axis2-&gt;SetBinLabel   (1-firstBin+kST,"");
	axis2-&gt;SetBinLabel   (1-firstBin+kZZ,"");
	axis2-&gt;SetBinLabel   (1-firstBin+kHWW,"");
	axis2-&gt;SetBinLabel   (1-firstBin+kNP,"");
	
	double xmin = h-&gt;GetXaxis()-&gt;GetXmin();
	double xmax = h-&gt;GetXaxis()-&gt;GetXmax();
	double ymin = h-&gt;GetYaxis()-&gt;GetXmin();
	double ymax = h-&gt;GetYaxis()-&gt;GetXmax();
	
	TGaxis* gaxis2 = new TGaxis(xmin,ymin,xmax,ymin,0,kNModes,511,"L");
	gaxis2-&gt;ImportAxisAttributes(axis2);
	
	TGraphAsymmErrors* cdf       = new TGraphAsymmErrors(kNModes);
	TGraphAsymmErrors* cdfprelim = new TGraphAsymmErrors(kNModes);
	TGraphAsymmErrors* d0        = new TGraphAsymmErrors(kNModes);
	TGraphAsymmErrors* d0prelim  = new TGraphAsymmErrors(kNModes);
	
	//  double msize = 1.25;
	double msize = 1.;
	
	cdf-&gt;SetMarkerColor(kBlue);
	cdf-&gt;SetLineColor  (kBlue);
	cdf-&gt;SetLineWidth  (2.0);
	cdf-&gt;SetMarkerStyle(21);
	cdf-&gt;SetMarkerSize (msize);
	
	cdfprelim-&gt;SetMarkerColor(kBlue);
	cdfprelim-&gt;SetLineColor  (kBlue);
	cdfprelim-&gt;SetLineWidth  (2.0);
	cdfprelim-&gt;SetMarkerStyle(25);
	cdfprelim-&gt;SetMarkerSize (msize);
	
	d0-&gt;SetMarkerColor(kBlack);
	d0-&gt;SetLineColor  (kBlack);
	d0-&gt;SetLineWidth  (2.0);
	d0-&gt;SetMarkerStyle(20);
	d0-&gt;SetMarkerSize (msize);
	
	d0prelim-&gt;SetMarkerColor(kBlack);
	d0prelim-&gt;SetLineColor  (kBlack);
	d0prelim-&gt;SetLineWidth  (2.0);
	d0prelim-&gt;SetMarkerStyle(24);
	d0prelim-&gt;SetMarkerSize (msize);
	
	double cdf_offset = 0.5; // 0.325;
	double d0_offset = 0.675;
	double xerr_width = 0.1;
	
	TBox* mybox = new TBox();
	mybox-&gt;SetFillColor(2);
	mybox-&gt;SetFillStyle(1001);
	
	TMarker* mymarker = new TMarker();
	mymarker-&gt;SetMarkerSize(msize);
	TArrow* myarrow = new TArrow();
	myarrow-&gt;SetLineWidth(2.0);
	
    for (int bin = firstBin; bin &lt; kNModes; bin++) {
		if (bin == kW || bin == kWg) {
			for (int i = 0; i &lt; 3; i++) {
				cdfvals   [bin][i] /= 0.10899;
				cdfvals_sm[bin][i] /= 0.10899;
				d0vals    [bin][i] /= 0.10899;
				d0vals_sm [bin][i] /= 0.10899;
			}
		}
		if (bin == kZ || bin == kZg) {
			for (int i = 0; i &lt; 3; i++) {
				cdfvals   [bin][i] /= 0.0366;
				cdfvals_sm[bin][i] /= 0.0366;
				d0vals    [bin][i] /= 0.0366;
				d0vals_sm [bin][i] /= 0.0366;
			}
		}
		
		if (cdfvals[bin][0]&gt;0) {
			if (0 &amp;&amp; log) {
				cout &lt;&lt;  bin &lt;&lt; " " &lt;&lt; cdfvals_sm[bin][0] &lt;&lt; " " &lt;&lt;  cdfvals_sm[bin][1]&lt;&lt; " " &lt;&lt; cdfvals_sm[bin][2] &lt;&lt; endl;
				//				mybox-&gt;DrawBox(bin+0.5-0.175, log10(cdfvals_sm[bin][0] -  cdfvals_sm[bin][2]),
				//							   bin+0.5+0.175, log10(cdfvals_sm[bin][0] +  cdfvals_sm[bin][1]));
				mybox-&gt;DrawBox(bin+0.175, log10(cdfvals_sm[bin][0] -  cdfvals_sm[bin][2]),
							   bin+1-0.175, log10(cdfvals_sm[bin][0] +  cdfvals_sm[bin][1]));
			} else {
// 				mybox-&gt;DrawBox(bin+0.15, cdfvals_sm[bin][0] - cdfvals_sm [bin][2],
// 							   bin+0.5,  cdfvals_sm[bin][0] + cdfvals_sm [bin][1]);
				mybox-&gt;DrawBox(bin+0.15, cdfvals_sm[bin][0] - cdfvals_sm [bin][2],
							   bin+1-0.15,  cdfvals_sm[bin][0] + cdfvals_sm [bin][1]);
			}
			
			if (cdfvals[bin][1] &gt; 0) {
				if (cdfisprelim[bin]) {
					cdfprelim-&gt;SetPoint(bin+1, bin+cdf_offset, cdfvals[bin][0]);
					cdfprelim-&gt;SetPointError(bin+1,xerr_width,xerr_width,  cdfvals[bin][2], cdfvals[bin][1]);
				} else {
					cdf-&gt;SetPoint(bin+1, bin+cdf_offset, cdfvals[bin][0]);
					cdf-&gt;SetPointError(bin+1,xerr_width,xerr_width, cdfvals [bin][2], cdfvals[bin][1]);
				}
			} else {
				double center = bin+cdf_offset;
				mymarker-&gt;SetMarkerColor(kBlue);
				mymarker-&gt;SetMarkerStyle(21);
				if (cdfisprelim[bin]) { mymarker-&gt;SetMarkerStyle(25); }
				mymarker-&gt;DrawMarker(center,cdfvals[bin][0]);
				myarrow-&gt;SetLineColor(kBlue);
				myarrow-&gt;DrawArrow(bin+cdf_offset, 0.84*cdfvals[bin][0], bin +cdf_offset, cdfvals[bin][0]*0.25,0.015, "&gt;");
			}
		}
		
// 		if (d0vals[bin][0]&gt;0) {
// 			if (log) {
// 				cout &lt;&lt;  bin &lt;&lt; " " &lt;&lt; d0vals_sm[bin][0] &lt;&lt; " " &lt;&lt; d0vals_sm [bin][1]&lt;&lt; " " &lt;&lt; d0vals_sm[bin][2] &lt;&lt; endl;
// 				mybox-&gt;DrawBox(bin+0.5, log10(d0vals_sm[bin][0] - d0vals_sm [bin][2]),
// 							   bin+0.85, log10(d0vals_sm[bin][0] + d0vals_sm [bin][1]));
// 			} else {
// 				mybox-&gt;DrawBox(bin+0.5, d0vals_sm[bin][0] - d0vals_sm[bin][2],
// 							   bin+0.85,  d0vals_sm[bin][0] + d0vals_sm[bin] [1]);
// 			}
			
// 			if (d0vals[bin][1] &gt; 0) {
// 				if (d0isprelim[bin]) {
// 					d0prelim-&gt;SetPoint(bin+1, bin+d0_offset, d0vals[bin][0]);
// 					d0prelim-&gt;SetPointError(bin+1,xerr_width,xerr_width,  d0vals[bin][2], d0vals[bin][1]);
// 				} else {
// 					d0-&gt;SetPoint(bin+1, bin+d0_offset, d0vals[bin][0]);
// 					d0-&gt;SetPointError(bin+1,xerr_width,xerr_width, d0vals[bin] [2], d0vals[bin][1]);
// 				}
// 			} else {
// 				double center = bin+d0_offset;
// 				mymarker-&gt;SetMarkerColor(kBlack);
// 				mymarker-&gt;SetMarkerStyle(20);
// 				if (d0isprelim[bin]) { mymarker-&gt;SetMarkerStyle(24); }
// 				mymarker-&gt;DrawMarker(center,d0vals[bin][0]);
// 				myarrow-&gt;SetLineColor(kBlack);
// 				myarrow-&gt;DrawArrow(bin+d0_offset, 0.84*d0vals[bin][0], bin +d0_offset, d0vals[bin][0]*0.25,0.015, "&gt;");
// 			}
// 		}
    }
	
	text = new TLatex(0,0, "");
	text-&gt;SetTextFont (62);
	text-&gt;SetTextAlign(21);
	text-&gt;SetNDC      ();
	text-&gt;SetTextSize (0.045);
	text-&gt;SetTextColor(kBlack);
	text-&gt;DrawLatex   (0.53,0.93,"CDF Run II, p#bar{p} at #sqrt{s} = 1.96 TeV");
	text-&gt;SetTextSize (0.025);
	text-&gt;DrawLatex   (0.89,0.0125,"M_{H}=160");
	
	cdf      -&gt;Draw("p");
	cdfprelim-&gt;Draw("p");
	//	d0       -&gt;Draw("p");
	//	d0prelim -&gt;Draw("p");
	
	// Setup legend
	//
	TH1F* sm = new TH1F("sm","dummy",1,0,1);
	sm-&gt;SetFillColor(2);
	sm-&gt;SetFillStyle(1001);
	sm-&gt;SetLineColor(152);
	
	TLegend* leg1 = new TLegend(0.63,0.65,0.75,0.88);
	leg1-&gt;SetBorderSize(0);
	leg1-&gt;SetFillStyle (0);
	leg1-&gt;SetTextSize  (0.040);
	leg1-&gt;AddEntry     (cdfprelim, "CDF Preliminary", "p");
	leg1-&gt;AddEntry     (cdf,       "CDF Published",   "p");
	//	leg1-&gt;AddEntry     (d0prelim,  "D0 Preliminary",  "p");
	//	leg1-&gt;AddEntry     (d0,        "D0 Published",    "p");
	leg1-&gt;AddEntry     (sm,        "Theory",          "F");
	leg1-&gt;Draw         ();
	
	c-&gt;SaveAs("XsecSumPlotFull.eps");
	c-&gt;SaveAs("XsecSumPlotFull.gif");
	
}
</PRE></BODY></HTML>
